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We consider an evolution operator for a discrete Langevin 
equation with a strongly hyperbolic classical dynamics and 
Gaussian noise. Using an integral representation of the evo- 
lution operator C we investigate the high order corrections to 
the trace of C n . The asymptotic behaviour is found to be con- 
trolled by sub-dominant saddle points previously neglected in 
the perturbative expansion. We show that a trace formula 
can be derived to describe the high order noise corrections. 



In the statistical theory of dynamical systems the de- 
velopment of the densities of particles is governed by 
a corresponding evolution operator. For a repeller, the 
leading eigenvalue of this operator C yields a physically 
measurable property of the dynamical system, the escape 
rate from the repeller. In the case of deterministic flows, 
the periodic orbit theory jl]]^] yields explicit and numer- 
ically efficient formulas for the spectrum of C as zeros of 
its spectral determinant (^]. 

On all dynamical evolutions in nature stochastic pro- 
cesses of various strength have an influence. In a series of 
papers (§-§ the effects of noise on measurable properties 
such as dynamical averages in classical chaotic dynami- 
cal systems were systematically accounted. The theory 
developed is closely related to the semi-classical h expan- 
sions [§-^o| based on Gutzwiller's formula for the trace 
in terms of classical periodic orbits jllj in that both are 
perturbative theories in the noise strength or Ti, derived 
from saddle-point expansions of a path integral contain- 
ing a dense set of unstable stationary points. The anal- 
ogy with quantum mechanics and field theory is made 
explicit in B where Feynman diagrams are used to find 
the lowest nontrivial noise corrections to the escape rate. 

In H we developed an explicit matrix representation 
of the stochastic evolution operator. The numerical im- 
plementation made it possible to reach up to order eight 
in expansion order, and the corrections to the escape rate 
were found to be a divergent series in the noise expansion 
parameter. This reflects that the corrections were calcu- 
lated (using the so called cumulant expansion) from other 
divergent quantities, the traces of the evolution operator 
C n {j. 

In [ 7) the focus was on the high order noise correc- 
tions for the special case of the first trace, Tr£. The 
asymptotics of the trace of the evolution operator were 



governed by sub-dominant saddles previously neglected 
in the expansion. 

In this paper we show that the high order noise cor- 
rections of TrC n are also dominated by sub-dominant 
saddles. These sub-dominant saddles can be treated as 
generalised periodic orbits of the system and we asso- 
ciate them with periodic orbits of corresponding discrete 
Newtonian equations of motion. Our key result is ( |40| ) 
where the high order noise corrections are converted into 
a trace formula. We give as a numerical example the 
quartic map considered in 

First we introduce the noisy repeller and its evolution 
operator. An individual trajectory in presence of additive 
noise is generated by iterating 



-1 = /On) + CT&, 



(1) 



where f[x) is a map, £ n a random variable with the nor- 
malised distribution and a parametrises the noise 
strength. In what follows we shall assume that the map- 
ping f(x) is one-dimensional and expanding, and that 
the £ ra are uncorrelated. A density of trajectories <f>(x) 
evolves with time on the average as 

</Wi(2/) = (£o <f> n ) (y) = J dx £(y,x)(f) n (x) (2) 

where the C evolution operator has the general form 

C(y,x) = 5 a (y - f(x)), 

S a (x) = j 'S(x- aO P (m = \v (f ) ■ (3) 

For the calculations in this paper Gaussian weak noise is 
assumed. In the perturbative limit, a — > 0, the evolution 
operator becomes 



C{x,y) 



1 



(y-/(»)r 



(4) 



The map considered here is the same as in our previous 
papers, a quartic map on the (0, 1) interval given by 



fix) = 20 



1 

16 



(5) 



Throughout the theory developed in previous works 

, the periodic orbits of the system played a major role. 

A periodic orbit of length n was defined simply by 



1 



Xj+l = f(xj), 

X n +1 = Xl- 



i = i, 



(6) 
(7) 



For a repeller the leading eigenvalue of the evolution 
operator yields a physically measurable property of the 
dynamical system, the escape rate from the repeller. In 
the case of deterministic flows, the periodic orbit theory 
yields explicit formulas for the spectrum of C as zeros of 
its spectral determinant ||. One of the most important 
goals of the theory related to stochastic evolution oper- 
ators is to explore the dependence of the eigenvalues v 
of C on the noise strength parameter a. The eigenvalues 
are determined by the eigenvalue condition 



F(a, v{a)) = det(l - C/v{a)) = 



(8) 



where F(a, 1/z) — det(l — zC) is the spectral determinant 
of the evolution operator C, which can be expressed as 



/ 7 

det(l - zC) = exp - 2] — Tr£ " 



(9) 



Equation (||) shows that noise dependence of the eigen- 
values of the evolution operator are very closely related 
to the noise dependence of the trace of C n , which shall 
be the object of study from now on. 
The trace of C n can be expressed as 



Tr£" 



(V27R7) 



1 / _ S 

dx\dx2.--dx n e ~^ , 



where 



1 ™ 



X n +1 = Xi. 



(10) 

(11) 

(12) 



In order to give a deeper insight on the forthcoming 
calculations we draw a correspondence between discrete 
Hamiltonian mechanics and our system, with the S de- 
fined above playing the role of the classical action. Ac- 
cording to (fijl) , the least action principle requires 



Xj - f(xj-i) - f'(xj)(x j+ i - f(xj)) = 0. (13) 
We define 



Pj ■= X J ~ f(x 3 -i), 



(14) 



the quantity corresponding to the momentum in the clas- 
sical mechanics. From (13) we obtain 



Xj+i = f{xj) +Pj+1, 



Pj+i 



f{x 3 y 



(15) 
(16) 



which are the equations corresponding to the classical 
Newtonian equations of motion. The generalised peri- 
odic orbits of length n are those orbits, which obey these 



equations and x n +i = x\, p n +i — Pn- Those generalised 
periodic orbits which have non-zero momentum will con- 
trol the asymptotic behaviour of the corrections to Tr£ n 
as we shall demonstrate later. The original periodic or- 
bits defined by (|^),(0) are those with zero momentum. 
The generalised periodic orbits with non-zero momentum 
and the original periodic orbits proliferate with growing 
n as suggested by Fig |l|. 
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FIG. 1. The sets of original and generalised periodic or- 
bits. Squares indicate original periodic orbits, dots indicate 
generalised periodic orbits, large symbols indicate orbits of 
length one, small symbols indicate orbits of length two. 

We introduce an integral representation of the noisy 
kernel, which will be of great use in the later calculations: 



C{x,y) 



1 («-/(x)) 2 
e 2^ 



2na 

27T 



(17) 



Using this new integral representation, 



Tr£ n 
1 

(2tt)" 



dk n dx n e~ 



or equivalcntly 



Tr£" : 
1 

(2tt)™ 



dk n I dp n J{p)e~ 



(18) 



. (19) 



where J(p) = D(x) / D(p) denotes the Jacobian. Since 

•f C n j}^lj = l k j P j 



we can reduce (|l?J) to 



(20) 



2 



Tr£" = 



r 2 n 

J dp n J(p)e^ A l[S( Pj ) 

~ 2a J(p) 



e 2 

where A denotes the Laplacian 
dp\ dpi 



(21) 



(22) 



Our object of study is the Taylor expansion of (^lj) in the 
noise parameter: 



Tr£" = ( Tr ' c ")jv cr 



2N 



N=0 

1 A 



jV 



(23) 
(24) 



vr- 



The iV-th power of the Laplacian in the equation above 
can be written as 



A-= £ 



. jll-jn } -dpl jl dp 



(25) 



where Sji is the Kronecker-delta. With the help of the 
multidimensional residue formula from complex calculus 



■•/(*) 



(2™) fe 

we obtain 

xS N yn 

'Z— /fe=l ■ 



(6 



(27ri)™2 



JV 



in.! +1 



E 



■(6 



(2ji)! 



\"fc+i : 



(2jn))! 



(26) 



h,-,3n=0 

J(p)dp 1 ...dp n 

2ji+i 2j„+r 

. Pi ---Pn 



(27) 



The contours are around the = points. The in- 
tegrals can be transformed back to contour integrals in 
the original Xj variables, and the contours will be placed 
around the original periodic orbits of the system defined 
by (§| — 0) , since it is these orbits which fulfil the pj = 
conditions. From now on we shall restrict our calcula- 
tions to the asymptotic large N limit. We will replace the 
summations in ( p7| ) by integrals and then use the saddle- 
point method to get a compact formula for (TrC n )N- We 
approximate the factorials via the Stirling-formula |l5| ] 
as 



(2j*)! ( 2 -ff k V^Tk 



Jfe! 



_ 2 2 Jk + l/2-jk e -J k 



(28) 



Using (|28|) and an integral representation of the delta 
function we get 



(Tr£ n ) 



x exp 



2- 



-JV 



N 



(2ni) n 2n 



E 



dt 



dx\...dx n 



n n n 

it(N - i*) + ( 2 ln 2 - !) E ^ + E ^ ln i k 



k=l 



k=l 



fc=l 



^ln(x fe -/(x fe _i))(2j fe + l) 



k=l 



(29) 



Now we replace jk with the new variables yk = and in 
the asymptotic (N large) limit approximate the summa- 
tions by yk with integrals by y k as 



(Tr/7% ~ 



(2m) n 2-K J 



dyi 



dy n / dt 
J j a 



dx\...dx n 



x exp 



it(N-Nj2yk) + N(21n2-l)J2vk 



fe=i 



fe=i 



+A^2/ fe ln(A^ fe ) + X>(*fc - f(x k ^))(2Ny k + 1) 

(30) 



fc=i 



fe=i 



We evaluate the ?/ integrals with the saddle point method 
to get 



(Tr£") 



2-iv+f 
w ~ (27r)fi"27r 

exp 





it '-' 






" e 2 



(31) 



Next we implement the saddle point method to the inte- 
gral in £ as well, asymptotically resulting in 



(Tr£ n ) 



N— 



(27V)! 



/v 



2 2JV +5(27r) 2 * 1 i™+ 1 M 
• i dx n e-( N+ ^ lo ^ s \ 



(32) 



The last step is to evaluate the contour integrals in the 
Xk variables. We deform the contours, until the saddle 
points are reached so the contours run along the routes 
of the steepest descent. The leading contribution comes 
from the saddle points, which fulfil the following equation 



[x* - /(z*_i) - (x* j+1 - f(x*))f'(x*)} = 0. (33) 



By comparing ( p3| ) and (|13|) one can see that the saddle- 
points are all generalised periodic orbits of the system. 
Since the contours ran originally around the orbits with 
zero momentum, these do not come into account as sad- 
dle points. The second derivative matrix is 



3 



(34) 



where D 2 S denotes the second derivative matrix of S 

d 2 S 



(D z S) t 



dxidxj 



(35) 



This would be the matrix to deal with if we would have 
taken the saddle point approximation of (^) directly. 
We reorganise the prefactor in (32) with the use of the 
Stirling formula [[13| and the result of the saddle point 
integration is written as 



(Tr£ n )jv * 



-N 



2m + ^&etD 2 S; 



(36) 



which is our main result. For n = 1 this formula gives 
back the result of as it should. 

Finally we draw the attention to the close connection 
between the generalised periodic orbits of the system and 
D 2 S. The stability matrix of a general periodic orbit is 
expressed as 



J, 



J — J\ ■ J% ■ J3... • J„ 

(f'{x k )- w ^ w f"{x k ) 



Pk 



f"{x k ) 



(37) 
(38) 



The determinant of D 2 S can be expressed with the help 
of the stability matrix as 



detD 2 S p = det( J p - 1) 
This way we reformulate ([56|) as 

r(AM 



(Tr£ n ) 



2n 



E 



-AflogSp 



(^V+f) f Vdet(l-J P )' 



(39) 



(40) 



where the summation runs over generalised periodic or- 
bits, with non-zero momentum. This is fully analogous 
to a trace formula and is our main result. 

Finally we turn towards testing our result obtained 
so far. In [Jt| we developed a contour integral method 
to calculate high order noise corrections to the trace of 
C. We showed that the agreement between the exact 
results and a formula which coincides with the ([36]) in 
the n = 1 case is very good. Now we step ahead and 
produce numerically high order noise corrections to the 
trace of C 2 . We shall start from ( ^jj) by transforming the 
integrals in p back to integrals in x as 



(Tr£ r ' 



1 



{2m) n 2 N 



E 

,...,j»=o 



(2ji)! (2 3n ))\ 



J 1 



x5 



dx\...dx n 



+1 ' 



(41) 



The contours at ( p7j ) were around the pj = points, so 
the contours above are placed around the original peri- 
odic orbits of the system, defined by @,(0). These con- 
tour integrals can be evaluated numerically. The Fig. ^ 
shows the ratio of (Tr£ 2 )^r obtained from (^6|) and eval- 
uated via the procedure described above as a function of 
N. 




FIG. 2. The ratio of (Tr£ 2 )iv calculated via the asymp- 



totic formula ( |36[ ) to its value computed by numerical inte- 
gration 

In summary we have studied the evolution operator 
for a discrete Langevin equation with a strongly hyper- 
bolic classical dynamics and a Gaussian noise distribu- 
tion. Using an integral representation of the evolution 
operator C we have revealed the asymptotic behaviour 
of the corrections to the trace of C n . This behaviour is 
governed by sub-dominant terms corresponding to terms 
previously neglected in the perturbative expansion, and a 
fully analogous trace formula can be derived for the late 
terms in the noise extension series of the trace of C n . 
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